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Abstract 

In recent years the theory of border coUision bifurcations has been developed for 

picccwisc smooth maps that are continuous across the border, and has been success- 
fully applied to explain nonsmooth bifurcation phenomena in physical systems. How- 
ever, many switching dynamical systems have been found to yield two-dimensional 
piecewise smooth maps that are discontinuous across the border. The theory for 
understanding the bifurcation phenomena in such systems is not available yet. In 
this paper we present the first approach to the problem of analysing and classify- 
ing the bifurcation phenomena in two-dimensional discontinuous maps, based on a 
piecewise linear approximation in the neighborhood of the border. We explain the 
bifurcations occurring in a physical system - the static VAR compensator used in 
electrical power systems, using the theory developed in this paper. This theory may 
be applied similarly to other systems that yield 2-D discontinuous maps. 



1 Introduction 

In recent years the bifurcations occurring in switching dynamical systems have 
been the subject of great interest, as it is known that many physical, en- 
gineering and biological systems actually embody continuous-time evolution 
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punctuated by discrete switching events. Under stroboscopic sampling these 
switching dynamical systems in general yield maps that have a discontinuity 
in the derivative along subspaces or "borderlines" that divide the phase space 
into two or more compartments. In these piecewise smooth (PWS) maps a 
new type of bifurcation, called border collision bifurcation (P E) |3l Sj) occurs 
when a fixed point collides with a borderline, resulting in a sudden change in 
the Jacobian matrix. Most recent studies on border collision bifurcations have 
been done on piecewise smooth maps that are continuous across the borderline 

dSJElElITD. 

It has been reported that many switching dynamical systems yield maps that 
not only have a discontinuity in the derivative, but also a discontinuity in the 
function itself. Discontinuity in the map arises if, in the n dimensional Poincare 
section there exists an (n — 1) dimensional manifold such that infinitesimally 
close points at the two sides of the manifold map to points that are far apart. 
A few electronic circuits exhibiting one- dimensional discontinuous maps was 
given in The bifurcation theory for such one-dimensional discontinuous 
maps have been developed (0 [10]), and have been successfully applied in ana- 
lyzing bifurcations in physical systems (II ip . However, in recent years there has 
been evidence that many physical systems yield two-dimensional (2D) discon- 
tinuous maps. For example, the classical impact oscillator yields a 2D discon- 
tinuous map if the impacting surface moves following a nonsmooth function 
(I12p . The static VAR controller used in electrical power systems (IT^ [Til) also 
has its dynamics given by a 2D discontinuous map. It has also been found that 
the dynamics of spiking bursting activities of real biological neurons ( ITS] [T6l) 
as well as some business cycle models fITTj) can be represented by 2D piecewise 
smooth discontinuous maps. Therefore, to explain the bifurcation phenomena 
in such systems, it is necessary to have a bifurcation theory for 2D piecewise 
smooth discontinuous maps. 

The bifurcation theory for ID and 2D continuous piecewise smooth maps 
is well developed ( !T8| [7]). The results related to the existence of period-1 and 
period-2 orbits in general n-dimensional PWS maps f5) have also been applied 
in practical systems. Bifurcation theory for ID discontinuous maps has been 
reported dHj). In the context of general n-dimensional discontinuous maps, some 
important results related to the existence of period-1 orbits (ITOll and period-2 
orbits f[T9l) have been published recently. For the special case of n = 2, these 
theories can be applied to obtain the conditions of existence of period-1 and 
period-2 fixed points. However in the context of a specific dynamical system, 
one is often interested in the asymptotically stable behavior, which cannot be 
inferred from the available theory. It has to be obtained through the analysis 
of the existence and stability of periodic orbits and their stable and unstable 
manifolds. Our investigations along this line showed that the two-dimensional 
discontinuous map is a source of incredibly rich dynamical behavior, created 
by complicated interactions between the stable and unstable manifolds of fixed 
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points and the line of discontinuity. In this paper we report some results of this 
investigation, and we hope that these will pave way for further investigation 
on the dynamics of two-dimensional discontinuous maps. 

The paper is organized as follows. In Section 2 we present an example of a 
physical system that yields a discontinuous map, and illustrate the peculiar 
dynamics exhibited by this system. As in the study of continuous piecewise 
smooth systems, we use the piecewise linear normal form representing the 
behavior of the system in the neighborhood of the border-crossing fixed point, 
which is introduced in Section 3. Then in Section 4 we develop a classification 
of the bifurcation scenarios observed in this map, depending on the type of 
the fixed points at the two sides of the border. 



2 A Physical Example: The Reactive Power Compensator 
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Fig. 1. The Static Var Compensator Circuit. The parameters are lo = 27r60 rad/s, 
Ls = 0.195 mH, i?, = 0.9 ml7, Lr = 1.66 mH, Rr = 31.3 ml7 and C = 1.5 mF. 

The system that we shall consider is the static VAR compensator used in elec- 
trical power systems, shown in Fig[Tj It is an inductor-capacitor combination 
(Lj. and C) in which the reactive power consumed by the inductor is controlled 
by back-to-back connected thyristor switches. The system is connected to a 
sinusoidal source u{t) at the input (representing the rest of the power sys- 
tem) through a transmission line of inductance Lg and resistance Rg. The two 
switches operate in the alternate half-cycles, and are turned on by applying 
pulses at a phase angle a with respect to the sinusoidal input u{t) = sin(co't). 
They turn off when the current through the respective switches reach a zero 
value. The system is described by the state vector x(t) of the current irit), 
capacitor voltage vdt) and the source current is(t). The dynamical equation 
during the ON state of either switch is 



x = Aonx + Bonm 



(1) 
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Fig. 2. The ir waveforms for two different initial conditions in Vc and is- 
where: 

^ -Rr/Lr l/Lr 

-l/C l/C 
-l/Ls -Rs/L 

During the OFF phase, v(^) is zero, and so the state vector is two dimensional, 
given by y = [vc is]'^- The corresponding dynamical equations are 

y = AoFFy + Boffm (3) 

where 

' l/C \ 10 
-l/Ls -Rs/Lj \1/L 



^OFF 



(4) 



We obtain the discrete map by observing the state vector stroboscopically, 
at each positive zero crossing of the input sinusoid. Since the stroboscopic 
observations are made when the switches are off, the discrete model must be 
two-dimensional. In this system, the switching logic imposes a discontinuity 
in the Poincare map. To illustrate, consider the waveforms of the inductor 
current v in Fig- [2| Suppose the switch is turned on at a phase angle a 
and the waveform of the current v is as shown. At a slightly different initial 
condition, the point of dip may reach a zero value, and so the instant of 
switch-off discontinuously changes as shown in the figure. If for a specific 
initial condition of Vc and ig the waveform just grazes the x-axis, then initial 
conditions at the two sides of this critical value map to widely separated points 
in the phase space. This implies that the resulting map is discontinuous. For 
a detailed derivation of the map, refer to (IT^t [T^ . 

Let us now examine the effect of this discontinuity on the system's bifurcation 
behavior. A typical bifurcation diagram with a as the variable parameter is 
presented in Fig. |3]^a). It is seen that for a > 90° and for a < 130°, the 
behavior is period- 1. But many peculiar dynamical transitions occur in the 
intervening parameter range. In such a system it is possible to locate the 
fixed point, irrespective of its stability, using a shooting method (12U|) . When 
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Fig. 3. (a) The bifurcation diagram of the static VAR compensator system, (b) 
Zoomed portion of the bifurcation diagram in the range (a G [98^^, 105'']) showing 
the occurrence of high-periodic orbits. 
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Fig. 4. The bifurcation diagram obtained by following the period-1 orbit using the 
shooting method. 



the bifurcation diagram is drawn by following the periodic orbit (Fig. |4j), it is 
found that there is no fixed point in those parameter ranges. Even though there 
is no fixed point, the orbit remains bounded and high periodic orbits occur in 
a specific order (see the zoomed portion in Fig. |3|^b)). The objective of this 
paper is to develop the theory necessary to explain such atypical bifurcation 
behavior in physical systems. 



3 The Normal Form 



Many hybrid dynamical systems can be represented in discrete-time by piece- 
wise smooth maps, given by equations of the form 



, q(x,y,u), (x,y) G Ra 
h{x,y,ij), {x,y) G Rb 



where /i is the bifurcation parameter and Ra and Rb are regions in the state 
space, divided by a borderline. We consider the class of maps that have the 



following properties: 



(1) The functions g and h are smooth (everywhere different iable), 

(2) The function / is discontinuous across the borderline, 

(3) The elements of the Jacobian matrix of / change discretely across the 
borderline, 

(4) The Jacobian elements are finite. 



Maps of the above properties have application in many physical and engineer- 
ing systems (the system described in Section [2] provides an example), and in 
the present paper we restrict our attention to such maps. 



The normal form for such a piecewise smooth system in the neighborhood of 
a fixed point on the border can be expressed as 
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where r^, tr are the traces and 5^, 6r are the determinants of the Jaco- 
bian matrices at the two sides of the border. The similar normal form of the 
piecewise smooth continuous map was derived in (P Ej). We have added the 
discontinuity / to obtain the normal form for the 2D discontinuous map. We 
shall call the two halves of the state space {x < 0} and {x > 0} as L and R 
respectively. 

In order to explain the bifurcations in a specific system, one has to obtain the 
eigenvalues of the fixed points at the two sides of the borderline, and from 
that, the values of the trace and the determinant. For example, in the case of 
the static VAR compensation system described in Section [2| the calculation 
of the Floquet exponents at a = 100° yielded tl = 1.4677, 6l = 0.6550, 
tr = 1.4677, and 6r = 0.6549. The explanation of the observed bifurcations 
have to be obtained by analyzing the behavior of the normal form map for 
these parameter values. 

The fixed points of the system ^ in both sides of the boundary are given by: 



L* 
R* 



fi + I -5R{fi + I) 



1 + 5r-tr l + 5_ 
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(7) 
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If the x-component of L* is negative, the fixed point exists. Else it does not. 
However, when the x-component of L* is positive, iterations from initial con- 
ditions in the left half are influenced by the "nonexistent" fixed point, which 
is called a "virtual" fixed point, and is denoted by L*. Similarly, when the 
x-component of R* is positive, the fixed point exists; else it is a virtual fixed 
point denoted by R*. 

The stability of L* and R* are determined by the eigenvalues 




We shall confine our studies to dissipative systems, so that \6l\ < 1 and 
\5r\ < 1. Moreover, we shall be concentrating on positive determinant sys- 
tems because most physical systems are observed to yield maps with positive 
determinant. For such a map, there can be four basic types of fixed points: 

(1) When -2y/S < T < 2^/6, both eigenvalues of the Jacobian are complex, 
with moduli less than 1, indicating that the fixed point is a spiral at- 
tractor. If r > (real part positive), it is a clockwise spiral, and if r < 
(real part negative), the motion is counter-clockwise. 

(2) When 6 < and —{1 + 6) < t < [l + S), both eigenvalues are real and 
are less than 1 in magnitude, causing the fixed point to be an attractor. 
If —(1-1-5) < r < —2^/6, then we have — 1 < A-t < 0, so that the attractor 
is a flip attractor (flips in both directions). If 2\/5 < r < (1-1-5), then 
< A+,_ < 1 and the fixed point is a regular attractor. 

(3) If r < —(1 -|- 6), then — 1 < A+ < and A_ < —1, so that the fixed point 
is a flip saddle (flips in both direction). 

(4) If r > (1-1-5), then A+ > 1 and < A_ < 1, so that the fixed point is a 
regular saddle. 

Earlier work of PWS maps (6) demonstrated a property of the normal form 
map that the unstable manifolds fold at every intersection with the x-axis, 
and the image of every fold point is a fold point. The stable manifolds fold at 
every intersection with the y-axis and the pre-image of every fold point is a 
fold point (jH [7]). In case of the discontinuous map, the same property holds; 
the only difference is that the folds are also discontinuous. 



4 Classification of discontinuous bifurcations 

We now classify the different types of discontinuous border collision bifurca- 
tions depending on the type of fixed points occurring at the two sides of the 
border. Unless otherwise stated, we shall study the bifurcations occurring in 
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the system (|6]) as the parameter /i is varied. Note that while in a continuous 
map (/ = 0) border colhsion occurs at a single value of fi, in the discontinuous 
case (/ 7^ 0) border collision occurs at two values n = —I and /i = 0. Note also 
that the structure of the map is different for positive and negative values of 
/. So in the subsequent sections we shall consider the bifurcations for positive 
and negative values of I separately. 

If a certain parameter combination exhibits a certain kind of bifurcation be- 
havior upon varying /i in one direction, then interchanging the parameter val- 
ues of the L and R sides will yield the same bifurcation behavior upon varying 
/i in the reverse direction, i.e., interchanging parameter values is qualitatively 
mirror-symmetric with respect to the bifurcation behavior. 

We now take up each case in the following subsections. 



4-1 Regular Saddle to Regular Saddle 
{tl>1 + Sl, tr>1 + 6r) 



Since both the fixed points are unstable, in most situations the orbits starting 
from all initial conditions go to infinity. However, when the discontinuity is 
negative, both L* and R* exist in the range /i G (0, — /). In such a situation, a 
trapping region can occur. 
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Fig. 5. (a) Structure of the stable and unstable manifolds in the range fj, S (0, —I) 
when L* and R* are regular saddles. The parameter values are = tr = 1.5, 
= = 0.3, / = —1, fi = 0.5. The trapping region is shown in white, (b) 
Bifurcation Diagram showing the existence of only a chaotic orbit in the trapping 
region. 

Let Ul and Sl be the unstable and stable manifolds of L* and Ur and Sr 
be the unstable and stable manifolds of R*, respectively. The discontinuously 
folded structure of the manifolds can be seen in Fig. ^a). Any initial state 
in L to the right of Sl diverges away from L* along Ul until it is mapped 
onto R. The state now comes under the influence of R* and is repelled back 
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along Ur. The state therefore gets locked between the two saddles, and this 
results in a bounded orbit. A similar behavior is seen for an initial state in R 
to the left of Sr. Since both sides are stretching, we observe a chaotic orbit. 
The chaotic attractor is stable as long as the unstable manifolds do not map 
to states outside the basin of attraction formed by 5*^ and Sr. This happens 
in the parameter range where 

-l{XuR-l) -1{1-Xsr) 

. _ ^-TR+Sn _ I-tr+Sr /n\ 

f^^ I-XSL , XuR-1 ' W 1~\sR I ^UL-l 

XuR is the eigenvalue of R* outside the unit circle, and Xsr is the eigenvalue 
of R* inside the unit circle. Xul and Xsl are the corresponding eigenvalues of 
L*. We will follow the same notations in the subsequent sections. Thus fii and 
fif are the parameter values at which boundary crisis occurs, i.e., the unstable 
manifolds Ur and Ul touch the stable manifolds Sl and Sr respectively. We 
see in Fig. |5]^a) that the unstable manifolds Ul and Ur meet the x-axis at A 
and B respectively while the stable manifolds Sl and Sr intersect it at C and 
D respectively. Therefore /Xj is the parameter value after which B lies towards 
the right of C and ;U/ is the value before which A lies towards the left of D, 
rendering the chaotic attractor stable. This range of occurrence of the chaotic 
attractor is seen in the bifurcation diagram of Fig. [sj^b). Note that if the map 
is continuous, i.e., if / = 0, the range of occurrence of the chaotic orbit goes 
to zero, satisfying the situation described in ([6]). 



4-2 Attractor (Regular/Flip/Spiral) to Regular Saddle 
(- (1 + 5l) < < (1 + 5l) , rR>l + 5r) 



The case of negative discontinuity: For /i < 0, both L* and R* are real and 
we have a period- 1 orbit. L* can be a regular/fiip/spiral attractor depending 
on the relationship between tl and 6l. The period-1 orbit is unique for < 
if L* is a regular attractor. However, high-periodic orbits as well as chaotic 
orbit may coexist with the period-1 orbit for values of /z slightly less than 
if L* is a spiral or flip attractor. When L* is very close to the y-axis, any 
trajectory approaching it has a possibility of going over to i? — which happens 
for spiralling or flipping orbits. 

For fi G (0,— Z), L* is a virtual attractor, and R* is a regular saddle. The 
stable manifold Sr of R* is responsible for forming the basin boundary, while 
the virtual attractor L* present in R draws the states in L within the basin 
of attraction towards R (see Fig. [6]^a), Fig. [7]^a), Fig. [Sj^a)). Any initial state 
in R to the left of Sr diverges asymptotically along the unstable manifold Ur 
until it is mapped onto L. It now comes back under the influence of L* and is 
drawn back to R. The trajectory therefore gets locked between the saddle and 



9 



the virtual attractor resulting in a bounded orbit. This orbit is stable before 
a parameter value Hc after which the unstable manifold Ur maps to points 
outside the basin boundary formed by Sr, i.e., a boundary crisis occurs at 
/i = /ic as /i is increased. The value of /ic depends on the slope of the unstable 
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Fig. 6. (a) The structure of stable and unstable manifolds when mupn > nisR and it 
is negative, tl = 1.0, tr = 1.0, Sl = Sr = 0.3, I = —1, fi = 0.2, and (b) Bifurcation 
diagram with tl = 1-0, tr = 1.5, 6l = 0.1, 6r = 0.2, / = —1. 

manifold after it first folds discontinuously at the This slope after the 

first fold is ruuRi = 5l/ i>^SR - tl) (Here \sr = Ar_). When m^/m > itisr, 
since both the slopes are negative, the two manifolds can never intersect (see 
Fig. Ha)). 



We see that Ur intersects the y-axis at A, and arbitrarily close points at the 
right and left of A map discontinuously to B and D respectively. The attractor 
becomes unstable due to boundary crisis when D moves to the right of C on 
the X-axis. This gives fic as 



(9) 



^ i2XsR -5r-1)1 

i^UR — ^Sr) 

which is positive. However, under the conditions depicted in Fig. [7](a) and 
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Fig. 7. (a) The structure of stable and unstable manifolds when muRi > msR but it 
is Positive, Tl = —1.0, tr = 1.0, Sl = Sr = 0.3, / = —1, fi = 0.2, and(b) Bifurcation 
diagram with ti = —1.0, tr = 1.5, 5l = Sr = 0.3, / = —1. 



Fig. |8[a), boundary crisis occurs when the segement oi Ur after the first fold 
intersects with Sr. We see from the corresponding Fig. [7](a) and Fig. |8]^a) that 
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Ur cuts the at B, and B maps to E. Thus the bounded orbit becomes 

unstable when E touches Sr. For this condition, /Xc is given by 



l^UR + (1 - ^ur) + tl^ur (1 - Xur) - \uRrB\ I 
[5r (1 + \ur) - 5l (1 - \ur) + t-lAc//? (1 - \ur) - ^ur'TrI 



(10) 
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Fig. 8. (a)The structure of stable and unstable manifolds when muRi < thsr, 
tl = 0.30, tr = 1.0, 5l = Sr = 0.3, / = — 1, /i = 0.2, and (b) Bifurcation diagram 
with TL = 0.30, TR = 1.5, 6l = 6r = 0.3, / = -1. 

As fi is increased through zero, the period-1 orbit vanishes through BCB and 
high periodic orbits of type L"i? come into existence (Fig. |6|^b) and Fig. ^h)) 
with period increment and Farey tree sequence. With further increase of /i 
symbol sequence gets reversed, i.e., the symbol sequence of type Li?" come 
into existence with period adding and Farey tree. But since tr > {1 + 6r), 
eventually the orbit becomes globally stretching, leading to chaotic dynamics 
at some value of /i. Depending on the value of tl a direct transition to chaos 
is also possible (Fig. [7|^b)). In both the cases chaotic attractor disappears due 
to a boundary crisis. 




Fig. 9. Bifurcation diagram with tl = 1.0, tr 



0.7, I 



-1. 



In the tl—tr parameter space, the domains of occurrence of each periodic orbit 
depends on the value of the determinants. For higher values of the determinant 
the regions overlap, resulting in coexistence of attractors and chaotic inclusions 
(see Fig. g. 
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For /i > — /, L* and R* are both virtual, which results in the entire state 
space being unstable because initial conditions in R diverge to infinity under 
the action of the unstable manifold (directed away from L to infinity) of the 
nonexistent fixed point R*, now located inside L, and initial conditions in L 
are also attracted to the side R due to the action of the nonexistent virtual 
attractor L*. 

For the case of positive discontinuity, we have a period-1 attractor in the case 
when both L* and R* are real (for < —I), else the entire space is unstable. 

The parameter space is symmetric upon variation of /i in the reverse direction 
when L* is a regular saddle and R* is an attractor, i.e., for tl > (1 + ^l) 
and — (1 + Sr) < tr < {1 + 6r). In that case a chaotic attractor is initially 
born out of a boundary crisis, slowly giving rise to increasingly stable orbits 
of non-monotonically decreasing periodicities, ultimately leading to a stable 
period-1 orbit. 



4-3 Regular Saddle to Flip Saddle 
{tl>1 + 5l, tr<-{1 + 6r)) 



The stable manifold Sl of the regular saddle in L is responsible for forming 
the boundary of the basin of attraction, while the fiip saddle in R is located 
inside this basin and is responsible for creating the attractor. This leads to 
many interesting phenomena as we shall see. 

The case of negative discontinuity: For yU < 0, when both saddles are virtual, 
we do not observe any attractor. For ^ > —I, L* is a regular saddle and R* is 
a flip saddle. All initial conditions in L converge onto the unstable manifold 
Ul- An initial condition to the right of the basin boundary formed by Sl 



converges onto the unstable manifold Ul and is mapped into R (see Fig. 10 [a)). 
The slope of the unstable manifold of R* is = —Sr/Xur = —Xr+ while 
the slope of the stable manifold is msR = —Xr-- The unstable manifold Ur 
folds discontinuously at the intersection with the x-axis and continues with 
a slope muR = SlXur/ {^r - tlXur) = ^lXr-/ {6r - tlXrJ). Since nisR > 
and muR < (as dictated by the range of parameters in question), the 
unstable manifold Ur must have a transverse homoclinic intersection with the 
stable manifold Sr after the first fold, which implies an infinity of homoclinic 
intersections and the existence of a chaotic orbit. For the regular saddle L*, 
we note that the unstable manifold Ul has a negative slope muL = —Xl+, and 
must, therefore, also have a heteroclinic intersection with the stable manifold 
Sr of the fiip saddle R*. The Lambda Lemma (l2Ti) implies that the unstable 
manifolds of the two sides will come arbitrarily close to each other, which 
makes the chaotic attractor robust 
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Fig. 10. (a) Homoclinic and heteroclinic intersections of the stable and unsta- 
ble manifolds of the flip saddle R* when tl = 1.5, = — 1.5, 5^ = 5_r = 0.3, 
/ = — 1, /i = 1.5. This figure and the subsequent ones showing the invariant 
manifolds have been generated using DsTool (|23| l25]l . (b) The chaotic at- 
tractor, and the basin of attraction created by the stable manifold of L* , for 
TL = 1.5, TR = -1.5, 5l = 5r = 0.3, I = -1, fi = 0.5. 



It is clear from the geometric structure that no point of the attractor can be to 
the right of point C. If C lies towards the left of D, the chaotic orbit is stable. 
If C falls outside the basin of attraction, it is an unstable chaotic orbit. From 
this we obtain the critical value of /i for which boundary crisis will occur: 



-Ac/. (1 





- 1) {5r - Xul^r) 



(11) 



which is the same as (10) taken in the reverse direction of /i variation. Note 
that for a continuous map the occurrence of the boundary crisis does not 
depend on /x, but in the discontinuous map it does. 



The same phenomenon occurs for < ^ < —I where R* is a virtual flip saddle, 
but its stable and unstable manifolds still exist in R and undergo homoclinic 
intersections resulting in chaos. The chaotic orbit occurs inside the basin of 
attraction formed by the stable manifold of L* (see Fig. [Io|(b)). 



The case of positive discontinuity: For /i > 0, L* is a regular saddle and R* 
is a flip saddle. In this case the same mechanism as discussed above causes a 



chaotic attractor, with the same stability condition as in (11). 



When there is a transition from a flip saddle to a regular saddle, i.e., when 
tl < — {1 + Sl) and tr > 1 + Sr, the bifurcation behavior is the same if /i is 
varied in the opposite direction. 
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4-4 Attractor to Attractor (Regular/ Flip /Spiral) 

- (1 + 6l) < < (1 + 5l) , - (1 + 5n) <tr<{1 + 6r) 



Once again, in this case, we shall be led to observe phenomena that do not 
occur in smooth maps or piecewise smooth continuous maps. 

The case of negative discontinuity: For both < and fi > —I, when one of 
the attractors is real and the other virtual, a stable period-1 orbit exists. If 
the fixed point is a spiral or flip attractor, high-periodic orbits may coexist 
with the period-1 orbit for values of n slightly less than zero or slightly greater 
than — /, because initial conditions close to the border may map to the other 
side before converging on the fixed point. 
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Fig. 11. (a) Bifurcation Diagram when both tl and r/j are positive, with tl = 1.0, 
6l = 0.2, tr = 1., = 0.2, and / = —1. (b) Bifurcation Diagram when tl is 
positive and tr is negative, with tl = 1.0, tr = —1.0, Sl = = 0.3, I = —1. 

For fi G (0, — /), both L* and R* are virtual attractors. Thus any initial state in 
L is drawn towards the virtual attractor L* situated in R. However as soon as it 
crosses the y-axis, it is drawn back towards the virtual attractor R* situated in 
L. The state therefore gets locked between the two virtual attractors and that 
results in a bounded orbit. The same behavior is exhibited by an initial state 
in R. This closed orbit may be high-periodic or chaotic. For < tl < {1 + Sl) 



and < tr < (1 + Sr), the bifurcation diagram in Fig. 11 'a) shows that as 
H is increased through zero period-1 orbit vanishes through border collision 
bifurcation and high periodic orbits of type L"i? come into existence with 
period inclusion and Farey tree sequence. With further increase of /i, the 
symbol sequence gets reversed, i.e., the symbol sequence of type LR^ comes 
into existence. Finally a period-1 orbit come into existence through BOB. 

We have seen that in a discontinuous map a parameter range can exist where 
both the fixed points are virtual, and if both the virtual fixed points are 
attracting in nature, high periodic orbits occur in specific sequence. 

However for higher values of the determinants {Sl,6r), the bifurcation struc- 
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Fig. 12. (a) Bifurcation Diagram when both tl and tr are positive, with 
tl = Tji = 1.4677, Sl = 6ji = 0.655, / = —1. (b) Bifurcation Diagram when 
both tl and tr are negative, with tl = tr = —1.0, 5l = = 0.3, I = —1. 

ture is more complex. Here we can observe period increment with coexistence 
of attractors. Notice that in Fig. [I2|^a) the parameters are the same as the 
ones obtained for the physical system considered in Section 2. The traces and 
determinants obtained from the physical system indicate that the local behav- 
ior of the system satisfies the condition being discussed in this section. The 



dynamical phenomena observed in Fig. ^ and in Fig. 12 a) are qualitatively 



the same, which explains the peculiar dynamical behaviors observed in the 
static VAR controller system. 
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Fig. 13. (a) Bifurcation diagram with respect to tr, for ti = 1.0, 6l = = 0.3, 
I = —1, /i = 0.5. (b) HomocUnic intersection of stable and unstable manifolds of a 
period-3 saddle. 

For < Tl < (1 + 6l) and — (1 + 5/?) < tr < 0, the high period orbits do 
not occur throughout the parameter range and a chaotic orbit is found to 
occur (see the bifurcation diagram in Fig. [ll|^b)) for high values of fi. This is a 
novel phenomenon — the occurrence of chaos in a system where both parts of 
the state space are contractive. The mechanism of the creation of the chaotic 
orbit becomes clear from the bifurcation diagram with th as the parameter 
(see Fig. [I3|(a)). As the parameter is reduced, the period-3 orbit's eigenvalues 
reach —1 and it undergoes a period doubling. The resulting period-6 orbit hits 
the border and turns into a six-piece chaotic orbit. The period-3 fixed point 
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is now a saddle, and Fig. 13 b) shows that its stable and unstable manifolds 
undergo a homoclinic intersection. This gives rise to the chaotic orbit. At a 
lower value of the parameter, a period-4 orbit starts to exist, initially coexisting 
with the chaotic orbit. At a specific parameter value, a boundary crisis occurs 
and the chaotic attractor disappears. At other parameter ranges, a similar 
phenomenon occurs for the other high-periodic fixed points. 

For — (1 + Sl) < tl < and — (1 + Sr) < tr < 0, i.e., when both the fixed 
points are flip attractors, a different situation occurs. In the parameter range 
fi G (0, — /) when no fixed point exists (both are virtual), only a period-2 orbit 
occurs (see Fig. [l2|b)). In fact, the range of occurrence of the period-2 orbit 
extends beyond this range, and the period-2 orbit coexists with the period-1 
orbit over a significant parameter range. Note that in the case of a continuous 
map, a stable period-2 orbit cannot occur when the map is globally contractive 

da ED- 

A symmetric behavior is exhibited for — (1 + 6r) < tr < and < < 
(1 + 6l) upon varying fi in the reverse direction. 
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Fig. 14. (a) Basin of attraction, for tl = 1.0, tr = 1.0, 5l = 5r = 0.3, fj, 
I = 1.0. (b) Bifurcation diagram with respect to fi, for tl = 1.0, tr - 
6l = Sr = 0.3, I = 1.0. 
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The case of positive discontinuity: In this case there is one real attracting fixed 
point in each of the parameter ranges /i < — / and /i > 0. For fj, e (0, — /) both 
L* and R* are attracting fixed points, and have their own basins of attraction. 
It is known that the basin boundary is normally formed by the stable manifold 
of a saddle fixed point. However, under the condition < < (1 -|- 6l) and 
< Tr < {1 + Sr) the basin boundary is formed by a completely different 
mechanism. In this case the segment of the borderline from —{fi + l) and — /i 
(the two preiterates of the origin) form a part of the basin boundary. The rest 
of the basin boundary is formed by successive preiterates of this segment, as 
can be seen in Fig. [I4|^a). This is due to the fact that in a discontinuous map, 
the borderline can act as a repellor. Moreover, its preiterates can also act as 
repellors since two points at the two sides of these line segments eventually 
map to points far apart. This mechanism gives rise to the bifurcation diagram 
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of Fig. [Tifb). 



4-5 Regular Attractor to Flip Saddle 

< Ti < (1 + , Tji<-{1 + 5r 



The Case of Negative Discontinuity: For < 0, L* is a regular attractor and 
R* is a virtual flip saddle, and the attracting fixed point enforces a stable 
period- 1 orbit since all initial conditions eventually lead to it. 

For /i > — /, L* is a virtual attractor, and hence lies in R. R* on the other 
hand is a flip saddle, and its stable and unstable manifolds undergo homoclinic 



intersections after the first fold (see Fig. 15 a)) to yield a chaotic attractor. 
After the first fold of Ur, its slope becomes 6lXur / {^r — tl^Ur) (a negative 
value), while the slope of Sr given by —\ur remains positive. All initial condi- 
tions in L are drawn to R due to the action of the virtual attractor L*, where 
they converge onto the chaotic attractor. This chaotic attractor extends into 
the parameter range < ;U < — /, when R* becomes virtual, but its stable 
and unstable manifolds still exist in R and undergo homoclinic intersections 
to yield a chaotic orbit. Hence, any initial condition in the entire phase space 
ultimately converges on the chaotic attractor. The chaotic attractor manifests 
for most of the parameter space, but for values of slightly greater than zero 
some high-periodic orbits are also stable. 
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Fig. 15. (a) Homoclinic intersections of stable and unstable man- 
ifolds of the flip saddle R* for a negative discontinuity, with 
tl = 1-0, tr = —1.0, = = 0.3, / = —1, ^ = 0.5. (b) Homoclinic in- 
tersections of stable and unstable manifolds of the flip saddle R* for a positive 
discontinuity, with tl = 1.2, tr = —1.5, 5l = Sr = 0.3, I = 1, fi = 0.5. 

The Case of Positive Discontinuity: For /i < — /, L* is a regular attractor and 
R* is a virtual flip saddle; for /i G (— /, 0), L* is a regular attractor and R* is a 
flip saddle; for /i > 0, L* is a virtual regular attractor and R* is a flip saddle. 
For the first two cases, the regular attractor in L causes a period-1 orbit to 
exist for all initial conditions. For the case when L* is a virtual attractor and 
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R* is a flip saddle, the stable and unstable manifolds of R* undergo homoclinic 
intersection (Fig. [I5|(b)) to yield a chaotic orbit. 



4-6 Spiral Attractor to Flip Saddle 
<tl< 2v^, tr<- 



R 



The Case of Negative Discontinuity: For /x < 0, L* is a spiral attractor and 
R* a virtual flip saddle, and the attractor causes a period-1 orbit. Whenever 
the spiral attractor is very close to the ?/-axis (which happens either when /i 
slightly less than zero or tl is negative), high-periodic orbits coexist with the 
period-1 orbit (see Fig. 16 a)). 
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Fig. 16. (a) Bifurcation Diagram with tl 
(b) Basin of Attraction with tl = 0.5, 
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High-periodic orbits continue to exist for /i G (0, — /) when L* is a virtual 
attractor and R* is a virtual flip saddle, because initial conditions in L are 
mapped to R due to the virtual attractor L* in R, and points in R are then 
mapped back to L due to action of the virtual flip saddle R*, thus leading to 
a bounded orbit, which manifests as a high-periodic or chaotic orbit (chaos 
occurs only when tl is positive). 

For fi > —I, L* is a virtual attractor while R* is a flip saddle. The virtual 
attractor L* maps all initial conditions in L to R, which subsequently flip to 
L. This may either lead to a stable high-periodic orbit or a chaotic orbit. The 
basin of attraction in Fig. [I6|(b) shows that there are high periodic attractors 
coexisting with a chaotic attractor when fi > —I. It is found that in much of 
the parameter range, the conditions for the occurrence of period-2 or period- 
3 are satisfled. The only difference between the cases when is positive 
(o < Tl < 2^/6L^ and when it is negative (^—2^/6L < tl < Oj (which basically 
implies a change in the sense of rotation of the spiral attractor L*) is that 
in the negative period-2 orbit is found to satisfy the conditions for 

stability and existence for the entire range of /i. 
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The Case of Positive Discontinuity: For /i < —I, L* is a spiral attractor while 
R* is a virtual flip saddle; for fi G (— /,0), L* is a spiral attractor and R* is a 
flip saddle; for > 0, L* is a virtual attractor while R* is a flip saddle. The 
first two cases will cause a stable period- 1 orbit to exist. For the last case, the 
virtual attractor L* maps all initial conditions in L to i?, which subsequently 
flip to L. This may either lead to a stable high-periodic orbit or a chaotic 
orbit. 



4.7 Flip Attractor to Flip Saddle 

(- (1 + Sl) <tl< -2v^, TR<-{1 + 5r)) 

The Case of Negative Discontinuity: For // < 0, L* is a flip attractor while 
R* is a virtual flip saddle; for /x G (0, — /), L* is a virtual flip attractor and 
R* is a virtual flip saddle; for ji > — Z, L* is a virtual flip attractor and R* 
is a flip saddle. For the first case a stable period-1 orbit exists. For the other 
two cases, all initial conditions in L flip to the other side of L* and land in 
R. Points in R diverge to infinity along the unstable manifold Ur of the fiip 
saddle (no matter whether it is real or virtual) , so that no attractor exists in 
the phase space. 

The Case of Positive Discontinuity: For /i < — /, L* is a flip attractor and 
R* is a virtual flip saddle. A stable period- 1 orbit exists but higher periodic 
orbit may coexist along with the period-1 orbit. For n G (— /,0), L* is a flip 
attractor and R* is a flip saddle, and a period-1 orbit exists and almost the 
entire space is stable. For > 0, L* is a virtual flip attractor while R* is a flip 
saddle, and no attractor exists. There is no basin of attraction because of the 
diverge-to-inflnity action of the unstable manifold of the flip saddle R* similar 
to the cases discussed earlier (since the attractor is virtual). 

The bifurcation behavior for the case tl < — {1 + Sl) and — (1 + Sr) < tr < 
—2\/5r is symmetric upon varying the parameter in the opposite direction. 



4-8 Flip Saddle to Flip Saddle (flipping along both directions) 
{tl < - (1 + 5i) , Tr < - (1 + 5lr)) 

No attractor exists in this case for any parameter value for both negative and 
positive discontinuity because of the action of the unstable manifolds of both 
the saddles — real or virtual. 
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5 Conclusion 



This paper contains the first attempt to analyze the bifurcation phenomena in 
two-dimensional discontinuous maps in terms of asymptotically stable behav- 
ior occurring for various parameter combinations. In this investigation we have 
used a piecewise affine approximation of the map in the neighborhood of the 
border, and then have partitioned the parameter space into regions depending 
on the types of the fixed points at the two sides of the border. For each case 
we have described the "typical" bifurcation phenomena in terms of the stable 
orbits. Where applicable, we have explained the occurrence and stability of 
chaotic orbits in terms of the structure of the stable and unstable manifolds 
of the fixed points. We do not claim to have done a complete study. Indeed, 
the dynamics of this system is so rich that it may take years to investigate all 
the possible situations. 

A few atypical features of the bifurcations in discontinuous maps have been 
found in this study: 

(1) The bifurcation behavior depends on the sign of the discontinuity; 

(2) There can be stable orbits even when no fixed point exists; 

(3) Chaos can occur even when each subsystem is contractive in nature; 

(4) Period-incrementing sequences are common. Each transition is "hard" in 
the sense that the state discontinuously jumps from one periodic orbit to 
another. 

(5) While in continuous maps, only stable manifolds of saddle fixed points 
can form a basin boundary, in discontinuous maps the line of discontinuity 
and its pre-images can also form a basin boundary. 

We have used the theoretical framework developed in this paper to explain 
the non-standard bifurcation behavior of the Static VAR Compensator. We 
have shown that the system yields a discontinuous map. As the firing angle 
is continuously varied, there comes a parameter range where no fixed point 
exists. In that situation, the state jumps discontinuously to high periodic orbits 
that occur in this range of the parameter, as per the prediction of our theory. 

We believe that this body of knowledge will help in understanding the non- 
standard bifurcations observed in a number of physical systems that yield 
discontinuous maps on discrete-time modeling. We also hope that this work 
will provoke further investigation on the dynamics of discontinuous maps. 
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